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An ideal atom laser would produce an atomic beam with 
highly stable flux and energy. In practice the stability is likely 
to be limited by technical noise and nonlinear dynamical ef- 
fects. We investigate the dynamics of an atom laser using a 
comprehensive one dimensional, mean-field numerical model. 
We fully model the output beam and experimentally impor- 
tant physics such as three-body recombination. We find that 
at high pump rates the latter plays a role in suppressing the 
high frequency dynamics, which would otherwise limit the 
stability of the output beam. 
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Optical lasers have had an enormous impact on sci- 
ence and technology, due to the intensity and coherence 
of the light they produce. We present a theoretical study 
of the output beam of the analogous matter wave device, 
known as the "atom laser" . These have been demon- 
strated in a number of laboratories, although so far only 
in an unpumped mode ■ 

Theoretical studies have highlighted both the similar- 
ities and the differences between optical and atom lasers 
[0-|l2|. The differences arise because atoms are more 
complex than photons: they have mass, giving a dif- 
ferent free space dispersion relation, and they interact 
with each other, producing strong nonlinearities. The 
latter generates complex dynamics, potentially compli- 
cating certain experimental measurements, of quantum 
noise for example. Nevertheless, the output beam dy- 
namics might be useful, in itself, as a probe of the exci- 
tations of the trapped condensate . 

Fundamentally, it is the coherence properties, such as 
linewidth and intensity correlations, that best capture 
the unique physics of lasers |l^Jl^. Coherence is also 
important for many practical applications, such as inter- 
ferometry. A classical nonlinear model of the laser as a 
noise driven van der Pol oscillator shows that the noise 
power, in both the phase and the amplitude, decreases 
in inverse proportion to the laser power [ p7[ . 

Schawlow and Townes showed that optical lasers' first 
order coherence, or linewidth, is ultimately limited by 
spontaneous emission, and also decreases in inverse pro- 
portion to the laser power [ p^ . Analogous results have 
been derived for atom lasers [|19"[. In practice, however, 
the Schawlow- Townes limit is not achieved, because the 
linewidth is limited by technical noise and by dynami- 
cal effects rather than by quantum noise. For example 
relaxation oscillations, due to the nonlinear interaction 



between the inversion and the light, are the primary de- 
terminant of the low frequency noise spectrum in many 
optical lasers . 

In this paper we focus on one aspect of the nonlinear 
dynamics of atom lasers: the frequency spectrum of the 
atom laser output beam density. In particular, we inves- 
tigate how the spectrum depends on the pumping rate. 
This spectrum describes the atom laser nonlinear dynam- 
ics. An experimentally measured spectrum would also 
include genuine "noise" due to physical processes that 
we do not model, such as fluctuations in the trapping 
potential. This investigation is motivated by the inverse 
relation between noise power and pumping rate in the 
optical laser. We choose to focus on the density because 
it can be measured experimentally using sensitive optical 
techniques such as phase modulation spectroscopy pl| ]. 

A conclusion of our work is that three-body recombi- 
nation plays a major role in the high frequency dynam- 
ics. As the pumping rate is increased the spectral peaks 
weaken and move to higher frequencies. This knowledge 
will facilitate atom laser applications, and measurements 
of their quantum noise. 




FIG. 1. Schematic representation of our atom laser model. 
The rings represent coils producing the trapping potential. 
The symbols are defined in the text. Atoms are injected into 
the non-condensed fraction. They can then either be lost or 
pump the condensate. They are coupled out of the condensate 
by a Raman transition and then fall in the gravitational field. 

Our atom laser model is schematically illustrated in 
Fig. 1. It consists of a dilute gas Bose- Einstein con- 
densate trapped by a highly anisotropic, effectively one- 
dimensional, potential. This is pumped from a reser- 
voir of non-condensed trapped atoms. The output atom 
laser beam is generated by a Raman transition to an 
untrapped atomic state. Our model builds on that of 
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Kneer et al. |ll|] by adding three-body recombination and 
a propagating output beam. In the following we will refer 
to the trapped Bose-Einstein condensate as the "conden- 
sate" , and to the output beam as the "beam" . We make 
the mean-field approximation in which the condensate 
and beam are described by the classical fields ^a{x,t) 
and '^b{x,t) respectively, which obey Gross-Pitaevskii 
type equations p^ . Our model is more complete than 
previous mean-field treatments [ pl] , p3| -[2^ since we in- 
clude pumping, output coupling, three-body recombina- 
tion, and explicitly model the dynamics of the pump 
reservoir, the condensate, and the beam. Although the 
extension to three dimensions is straightforward, we do 
not report it here. To describe the pumping we use a 
phenomenological model, which mimics the pumping of 
an optical laser [|ll|. It only depends on the total num- 
ber of atoms in the uncondensed component, Nu{t), not 
on its spatial structure. The condensate atoms are cou- 
pled, by a reversible Raman transition [p| j2^j27|] , to an 
untrapped electronic state forming the atom laser beam. 
The Raman transition imparts a momentum kick hk to 
the out-coupled atoms. The beam evolves under the in- 
fluences of gravity and atom-atom interactions, which are 
dominated by those due to the condensate. We have only 
considered the particular case for which the Raman cou- 
pling is tuned to transfer atoms starting from the center 
of the harmonic trap | p3[ |. Hence the output beam over- 
laps with the lower half of the condensate. Experimen- 
tally, three-body recombination is well described as a loss 
process occurring at a rate proportional to the cube of the 
local density We follow Kagan et al. in incor- 
porating three-body recombination into the atom laser 
Gross-Pitaevskii equations. 

Our model is defined by the following dimensionless 
equations: 



dt 



r - juNu - ipNuNa, 



2 dx"^ 2 



(1) 



dt 



2 



Uab^b\^a\'' -llr-^b{\^at+\'^b\ 
^ a ■ 



The model is made dimensionless using the character- 
istic trap length / = {H/ujmY^'^ and angular frequency 
uj, with m the atomic mass. Hence time i, position x, 
and the fields are measured in units of w^^, I and 
respectively. Experimentally reasonable values of the pa- 
rameters are: a trap frequency u; « 125 Hz pO|, and the 



atomic mass of sodium m = 3.8 x lO"^^ kg These give 
a time scale of w""'^ = 8 ms, and a length scale of / = 4.7 
fim. 

Atoms are injected into the uncondensed fraction at 
the rate r, are lost at the rate 7ii-/V„, and pump the con- 
densate at the rate jpN^Na, where Na = j^apc'a; 
is the total condensate population. Ua and Ub are the 
intra- and Uab the inter- species two-body interaction co- 
efficients. The dimensionally correct coefficients can be 
written as Ua,b,ab = 47ras/(?i/wm)^^2 with as the appro- 
priate s-wave scattering lengths . We have previously 
found that three-body recombination, with coefficient 7^, 
is necessary for the system to reach a quasi-stationary- 
state on a timescale comparable to experimental conden- 
sate lifetimes The three-body recombination rate 
7r = 10"'' is chosen to produce Na ~ 3500 — 9500 con- 
densed atoms at steady state, as the pump rate varies 
from r = 200 to 1600. For simplicity we assume that the 
inter- and intra-species three-body recombination rates 
are the same. G is the scaled gravitational acceleration 
Gravity has not been explicitly included for the 
trapped atoms, since it is equivalent to a spatial shift of 
the trap potential. The Raman coefficient is a func- 
tion of the Raman laser amplitude and detuning |27|. 

Our numerical method is a split-step Fourier method 
psf with a fourth order Runge-Kutta in place of the usual 
nonlinear step. We have absorbing boundaries at the 
ends of the spatial grids to prevent reflections. 
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FIG. 2. Densities of the condensate (dashed) and beam 
(solid) versus position. All quantities are dimensionless. The 
plotted beam density is 100 times the true value. These re- 
sults are a numerical solution of Eqs. (jl|) at t = 200 after 
starting from a seed condensate. Parameters are: r = 1600, 
Ua = Ub= 0.02, Uab = 0.01, G = 12, 7H = 0.5, 7« = 7p = 0.1, 
7r ~ 10~^ and k = 5. 

Typical results for quasi-stationary-state spatial pro- 
files of the condensate and beam densities are shown 
in Fig. 2. For our parameters, the condensate shape is 
well described in the Thomas-Fermi approximation. Note 
that the plotted output beam density has been scaled up 
by a factor of 100. The flux in the beam is about 90 
atoms per unit dimensionless time, or about 1.1 x 10^ 
atoms per second with w = 125 Hz. The prominent spa- 
tial oscillations of the output beam, occurring within the 
extent of the condensate, have been observed in previ- 
ous work and are related to the oscillatory Airy 
type eigenf unctions of particles in the linear gravitational 
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FIG. 3. The output beam density, |^'5(xo,t)P, as a func- 
tion of time at position, xq ~ —45, for different pump rates: 
from bottom to top r = 200, 400, 800, and 1600. All quanti- 
ties are dimensionless. Other parameters are as for Fig. |^. 

In Fig. ||we show the output beam density |^'fc(xo,t)P 
as a function of time at a fixed position far from the 
trapped condensate. The time span has been chosen 
to be longer than that of the transient dynamics asso- 
ciated with the growth of the condensate from a seed, 
and short enough to be experimentally accessible. Quasi- 
stationary, nearly periodic dynamics develop, becoming 
approximately harmonic for large pump rates r. We have 
verified that the quasi-stationary dynamics is indepen- 
dent of the particular form of the seed used to initiate the 
condensate growth. The fundamental oscillation of the 
beam, at higher pump rates, arises from the well known 
Kohn mode ||2^ , corresponding to a rigid "sloshing" back 
and forth of the condensate at the trap frequency. 



100 
50 





r.400 


il ...I.. 





■e "o 10 

<100| ^ 



I 50 
o 

a. „ 



Angular Frequency 
FIG. 4. Frequency spectra of the time dependent density 
of the fields in Fig. pi All quantities are dimensionless. The 
spectra are those of the time series from t = 112.5 to t = 175. 
The pump rates from top to bottom are r — 200, 400, 800, 
and 1600. The normalization has been chosen so that the 
largest peak height (that near 20 frequency units for r=200) 
is 100 units. The peaks just below 50 frequency units, for 
r = 1600, are almost invisible on this scale. 

Fig. ^ shows the frequency spectra of the density time 



series in Fig. |3[ and is the central result of this paper. The 
spectra are those of the final, quasi-steady, parts of the 
time series. The frequencies are well approximated by the 
eigen-frequencies of an uncoupled ID condensate. Kneer 
et al. [|ri| found these to be = uj^/n{n + l)/2 where 
n — 1,2,3,..., with the corresponding spatial modes 
given by the Legendre polynomials P„(x). The Kohn 
sloshing mode at w is n = 1, and n = 2 is the breathing 
mode at V^uj. In this mode the width of the condensate 
oscillates, with a corresponding density oscillation. 

Fig. ^ shows that as the pumping rate r increases, the 
slow beam oscillations change from the frequency of the 
condensate breathing mode to that of the sloshing mode. 
This is because the condensate itself changes from breath- 
ing to sloshing at higher pumping rates. As Kneer et al. 
|lT| noted, the sloshing mode is spatially asymmetric and 
hence must be excited by a spatially asymmetric pertur- 
bation, such as the output beam. This is consistent with 
the observed increase of the sloshing mode power with 
the density of the output beam, that is with pumping 
rate. On the other hand, mode damping due to three- 
body recombination increases rapidly with density. We 
have verified that in a simple model of the Kneer et al. 
im form, the breathing mode decays with pumping rate. 

Although the frequencies in the beam are primarily 
determined by the condensate, there is a complex rela- 
tionship between the power spectra of the condensate and 
of the beam. In particular, the beam spectrum reflects 
the spatial dependence of the condensate spectrum. Con- 
sider, for example, the condensate Kohn sloshing mode. 
The magnitude of the local density change varies with 
the spatial derivative of the density profile. It is maxi- 
mum where the condensate has greatest slope, and zero 
at the center where the slope is zero. That the Kohn 
mode frequency is prominent in the power spectrum of 
the output beam indicates that the coupling between the 
condensate and the beam is not confined to the center of 
the condensate. Similar observations apply for the high 
frequency components of the power spectrum, the en- 
velopes of which differ in detail between the condensate 
and the beam. 

As the pumping rate increases. Fig. 4 shows that the 
frequency of the high frequency group increases and its 
power decreases rather dramatically. A decrease of spec- 
tral power with pumping rate, due to gain saturation, is 
characteristic of optical lasers. Since three-body recombi- 
nation contributes to gain saturation in atom lasers, we 
have investigated its role in a simplified pumped atom 
laser model similar to that of Kneer et al. [|ll|. That 
is, with a spatially uniform phenomenological loss rather 
than an explicitly modelled output beam. Without three- 
body recombination the condensate density power spec- 
trum simply rolls off at high frequencies. Adding three- 
body recombination produces a high frequency spectral 
group, such as seen in Fig. 4. Hence it is an important 
factor in the dynamics of the system. 
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In conclusion, we have explored some of the rich and 
complex behavior of the atom laser with a model incor- 
porating important experimentally relevant physics, such 
as three-body recombination. The latter was found to 
strongly influence the spectrum of the atom laser output 
beam. This work is a step towards the future goal of un- 
derstanding and measuring the quantum noise properties 
of the atom laser. 
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